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ABSTRACT 

We present results from two Chandra imaging observations of Ross 154, a nearby flaring M dwarf 
star. During a 61-ks ACIS-S exposure, a very large flare occurred (the equivalent of a solar X3400 
event, with Lx = 1-8 x 10 30 ergs s _1 ) in which the count rate increased by a factor of over 100. The 
early phase of the flare shows evidence for the Ncupcrt effect, followed by a further rise and then a 
two-component exponential decay. A large flare was also observed at the end of a later 48-ks HRC-I 
observation. Emission from the non-flaring phases of both observations was analyzed for evidence 
of low level flaring. From these temporal studies we find that microflaring probably accounts for 
most of the 'quiescent' emission, and that, unlike for the Sun and the handful of other stars that 
have been studied, the distribution of flare intensities does not appear to follow a power-law with 
a single index. Analysis of the ACIS spectra, which was complicated by exclusion of the heavily 
piled-up source core, suggests that the quiescent Ne/O abundance ratio is enhanced by a factor of 
~ 2.5 compared to the commonly adopted solar abundance ratio, and that the Ne/O ratio and overall 
coronal metallicity during the flare appear to be enhanced relative to quiescent abundances. Based on 
the temperatures and emission measures derived from the spectral fits, we estimate the length scales 
and plasma densities in the flaring volume and also track the evolution of the flare in color-intensity 
space. Lastly, we searched for a stellar-wind charge-exchange X-ray halo around the star but without 
success; because of the relationship between mass-loss rate and the halo surface brightness, not even 
an upper limit on the stellar mass-loss rate can be determined. 

Subject headings: stars: coronae — stars: individual (Ross 154) — stars: late-type — stars: mass loss — 
X-rays: stars 



1. INTRODUCTION 

Apart from the Sun, the best opportunities to study 
stellar magnetic phenomena at extremely low levels are 
provided by nearby stars. Of particular interest are the 
M dwarfs, whose lower photospheric temperatures and 
often higher coronal temperatures represent significantly 
different conditions from those of the solar atmosphere. 

M dwarfs tend to be more active than F-G dwarfs, 
with approximately half of the former (the dMe stars) 
showing emission in Ha, the hallmark of flaring activ- 
ity. Despite having surface areas only several percent as 
large as the Sun's, dMe stars have X-ray emission that 
is usually comparable to or larger than the solar X-ray 
luminosity; the flaring coronae of M dwarfs would thus 
appear to present a much more dynamic environment 
than offered to us by the Sun. These conditions provide 
a means for determining how stellar activity differs from 
the solar analogy, with the ultimate goal of underpinning 
observational similarities and differences with the funda- 
mental physics needed to describe the various activity 
phenomena on display. 

Here we present analyses of two Chandra X-ray Obser- 
vatory observations of the M3.5 dwarf Ross 154, which 
counting Proxima Cen as pa rt of the a Cen system 
Wertheimer fc Laughlinl l2006t )) is the seventh nearest 
stellar system to the Sun. The original motivation for our 
primary observation, with the ACIS-S detector, was to 
constrain the stellar mass-loss rate by searching for X-ray 
emission from charge-exchange collisions of its ionized 
wind with the surrounding interstellar medium (ISM); 
Ross 154 is one of the few stars with a combination of dis- 



tance and mass-loss rate that might be amenable to this 
technique using current observing capabilities (sec SjSJ). 
During the observation, the star underwent an enormous 
flare during which the X-ray photon count rate rose two 
orders of magnitude above the quiescent value. These 
X-ray data offer a rare glimpse of the time evolution of 
flaring plasma while simultaneous optical monitoring ob- 
servations using the Chandra Aspect Camera Assembly 
provide insights into the response of the chromosphere 
and photosphere during the event. A second observation 
using the HRC-I detector provides only temporal infor- 
mation, but with superior statistical quality. 

In Sj2] wc summarize what is known about Ross 154. 
Sections [3] and |4] describe the X-ray and optical obser- 
vations and data reduction. The analysis of the ACIS 
X-ray spectra is described in $5l and temporal analyses 
of the photon event lists and large flare light curve are 
presented in §§\ and SJ3 respectively. Finally, we discuss 
our attempts to deduce the mass- loss rate of Ross 154 in 
^8] before summarizing our findings in i}9l 

2. THE TARGET 

Ross 154 (Gleise 729, V 1216 Sgr) i s a flaring M3.5 
dwar f and lies at a distance of 2.97 pc (jPerrvman et alJ 
119971) towards the Galactic Center (I = 11.31°, b = 
-10.28°; RA = 18:49:49, Dec. = -23:50:10). Its 
X-ray luminosity, estimated at 6.0 x 10 27 ergs s" 1 
based on Rontqen Satellite (ROSAT) measurements 
(|Hunsch et al.1 119991 ). is modest for active M dwarfs, 
which range up to nearly 1 30 ergs s -1 in quiescence. It s 
\og(L x /L hol ) value is -3.5 (|Johns-Krull fc Valenti|[l996D . 
somewhat below the saturation limit of ~ —3.0 for ac- 
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tive late-type stars ()Agrawal. Rao, fc Sreekantanl 119861 : 
iFleming et al.lfl993l ). RS CVn stars, which are close bi- 
nary systems — often tidally locked — with high rotation 
rotation rates and flaring activity, have lum inosities up 
to - 10 32 ergs s- 1 (|Kellett fc Tsikoudilll997[ ). 

Ross 154 is a moderately fast rotator, with v sinz = 
3.5 ± 0.5 km s _1 measured by iJohns-Krull fc Valentil 
(|l996t ). indicating an age of less than 1 Gyr. Those au- 
thors also estimate a magnetic field strength of 2.6 ± 0.3 
kG with a filling factor / = 50 ± 13%, somewhat 
weaker than the ~4-kG fields with / > 50% measured 
for three other more X-ray-luminous and faster-rotating 
dMe stars: AD Leo (Saar & Linsky 1985), EV Lac 
(Johns-Krull & Valenti 1996), and AU Mic (Saar 1994; all 
three stars). Its photosphcric mctallicity is roughly half- 
solar based on the estimated iro n abundance of [Fc/H]^ 
-0.25 reported bv lEggenl ((1996), and the presence of op- 
tical emission lines from Fe I, Si I . and Ca I indicates 
a sur prisingly cool chromosphere (jWallerstein fc Tvagil 
12001 . 

Ross 154 has been detected by a number of high-energy 
observatories — Einstein (Agrawal et al. 1986; Schmitt et 
al. 1990), ROSAT (Wood et al. 1994; Hiinsch et al. 1999), 
and the Extreme Ultraviolet Explorer (EUVE; Bowyer et 
al. 1996; Lampton et al. 1997; Christian et al. 1999) - 
but was never studied in any detail, and no significant 
flares were seen in the relatively short exposures obtained 
by those missions. 

3. THE OBSERVATIONS 

Two sets of Chandra data were analyzed. The first ob- 
servation was conducted for 60625 s beginning on 2002 
September 9 at 00:01:20 UT {Chandra time 147916880). 
The primary objective was to look for a faint stellar-wind 
halo around the star so no grating was used, despite the 
brightness of the coronal emission and the expectation of 
severe pileup in the source core. This measurement (Ob- 
sld 2561) used the ACIS-S3 backside-illuminated CCD 
chip in Very Faint (VF) mode with a subarray of 206 
rows that allowed a short 0.6-s frame time to be used, 
thus reducing the degree of event pileup. 

Standard Chandra X-ray Center pipeline products 
were reprocessed to level 2 using the Chandra Interac- 
tive Analysis of Observations (CIAO 1 ) software version 
3.4 with calibration data from CALDB 3.3.0.1, which 
applies corrections for charge-transfer inefficiency in all 
ACIS CCDs and for contamination build-up on the ACIS 
array. Eight weak secondary sources were then removed 
from the source field (Fig. [1^,) and VF-mode filtering 2 
was applied to reduce the background, except near the 
core of the primary source as explained in 2J No back- 
ground flares were observed in the data. 

Examination of the X-ray light curve (Fig. ^ after ex- 
cluding the piled-up core of the target reveals three tem- 
poral phases in our observation: a quiescent phase; a 
very large flare; and the flare decay. We also obtained 
a simultaneous optical light curve using data from the 
Chandra aspect camera. 

The second observation (Obsld 8356) occurred on 
2007 May 28 beginning at 4:19:18 UT (Chandra time 



1 |http : //exc .h arvard. ed u/ciaoT] 

2 See ~. Vikhlinin (2001 
http : //exc .harvard. edu/ciao/threads/aciscleanvf / 
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Fig. 1. — ACIS data extraction regions. Panel (a) shows the 
annulus used for Quiescent spectra data around the primary source, 
small circles around weak secondary sources, and two rectangular 
background regions; also note the bright readout streak running 
approximately horizontally through the primary source. Panel (b) 
shows the spectral extraction regions for the Quiescent, Flare, and 
Decay phases; the spatial scale is twice that of panel (a). Spectral 
extractions exclude the piled-up core and the readout streak, the 
latter of which may have a slightly different detector gain from 
non-streak regions. Streak data, however, are included in temporal 
analyses. 
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Fig. 2. — Optical and X-ray light curves from the ACIS ob- 
servation, with time bins of 20.5 and 100 s, respectively (500 s in 
X-ray detail). Zero time corresponds to the beginning of the ob- 
servation at 2002 September 9 at 00:01:20 UT. X-ray light curve 
is for E = 250 — 6000 eV with a spatial filter comprising an an- 
nulus of radii 6 and 60 pixels (to exclude the piled-up core) plus a 
6-pixel-wide box along the readout streak. The annulus used for 
the Quiescent Detail light curve had radii of 4 and 40 pixels. In 
the absence of pileup, so that core events could also be used, the 
counting rate would be ~25 times higher. Note the small optical 
flare at the very beginning of the X-ray flare. 



296713158) and ran for 48511 s, excluding two 32.8-s 
periods that did not meet Good Time Interval criteria. 
This HRC-I calibration measurement was made in paral- 
lel with a measurement of the ACIS background, in which 
ACIS is placed in a "stowed" position where it is both 
shielded from the sky and removed from the radioactive 
calibration source in its normal off-duty position. With 



Ross 154 X-Ray Flares 



3 




Fig. 3. — Image of the HRC-I off-axis observation. The inner 
ellipse marks the source region used for microflaring analysis. The 
background region used to construct the light curve in Fig. [8] lies 
outside the larger ellipse. 

ACIS in the stowed position, the HRC-I is 75 mm away 
from its nominal on-axis position but can be used for off- 
axis observations to calibrate the Chandra point spread 
function (PSF) and measure small-scale gain uniformity. 
When ACIS is collecting data in its standard telemetry 
mode the HRC can operate simultaneously, albeit with 
severe telemetry restrictions, in its Next-In-Line (NIL) 
mode with a limit of ~3.5 counts s _1 . Ross 154 is one 
of a handful of isolated sources with a counting rate and 
other parameters suitable for such NIL-mode calibration 
observations. 

The HRC-I observation was made 25.62' off axis (Y 
offset = -11.83', Z offset = 22.73') using a 10-tapxlO- 
tap (5.61' x 5.61' ) window of the detector, large enough 
to encompass the out-of-focus source image and also pro- 
vide a suitably large region to determine the background 
level (see Fig. Within the windowed detector re- 
gion the background accounted for ^2.6 counts s _1 and 
the total rate during source quiescence was ~3.0 counts 
s _1 . At least one large flare exceeded the telemetry limit. 
The HRC-I has essentially no spectral information and 
so these data were used only for the microflaring study 
presented in SJBJ Until then our discussion will refer ex- 
clusively to the ACIS observation. 

4. ACIS DATA EXTRACTION 

4.1. Optical Monitor Data 

Although the capability is rarely used, one or more 
of the eight Chandra Aspect Camera Assembly (ACA) 
"image slots" can be assigned to monitor selected sky 
locations, with a practical faint limit of roughly V = 12. 
For comparison, the guide stars used for aspect determi- 
nation and to correct for spacecraft dither generally have 
magnitudes between V — 6 and 10, and Ross 154 has a 
magnitude of V = 10.95. The faint limit is largely deter- 
mined by pixcl-to-pixel changes in the optical CCD dark 
current, with temporal variations on scales of typically a 
few ks. 3 The effects of these background variations are 
significant in our observation, but were largely removed 
by averaging data collected when the source centroid 
dithered more than ~3.5 pixels away from the pixel un- 

3 Chandra Proposers' Observatory Guide, section 5.8.3, at 
http : //exc . harvard . edu/proposer/POG/ 



der consideration and then subtracting the average bias 
within each pixel. Quantum efficiency nonuniformity was 
mapped by analyzing data from times when the source 
centroid was within 0.1 pixels of the center of the peak 
pixel; the flat field corrections thus derived were never 
more than 3%. 

In addition to having a low signal level, the optical 
monitoring of Ross 154 suffered a loss of tracking halfway 
through the quiescent phase of the observation, causing 
the star to wander within its 8 x 8-pixel ACA window 
(with integration time of 4.1 s per frame) as the space- 
craft dithered with a period of 707.1 s in pitch and 1000.0 
s in yaw, causing semi-periodic dips of up to 14% in the 
light curve as some of the source flux was lost near the 
edges of the window. Various flux correction schemes 
were tried, but the best results were obtained by sum- 
ming the background-subtracted counts within a 5 x 5 
box centered on the highest-count pixel and discarding 
all data from cases where the 5x5 box did not fit within 
the 8x8 ACA window. The resulting optical light curve 
is shown in the top panel of Figure^ revealing a 0.15- 
magnitude flare and a smaller precursor; these are dis- 
cussed in more detail in the context of the X-ray flare 
behavior in [J7] 

4.2. X-Ray Data 

After removing secondary sources from the X-ray data, 
a background region was obtained by excluding the 
dithered edges of the field and a 300-pixel-wide (147") 
region around the source (Fig. [IJl). Even with the large 
exclusion region around the main source, the flare was so 
intense that excess counts could still be seen in the back- 
ground light curve so we excluded 2000 s around the flare 
peak. The remaining background data show no evidence 
for temporal variability. 

Ross 154 is such a bright X-ray source that even dur- 
ing quiescence it produced a strong CCD readout streak 
and the core was very heavily piled-up. To obtain undis- 
torted spectra we excluded the core, with a different ra- 
dius for the quiescent, flare, and decay phases of our 
observation (see Table [1} , based on studies of spectral 
hardness ratios and VF-filtering losses as a function of 
radius. The hardness ratio method is based on the fact 
that piled-up spectra shift lower energy flux to higher 
energies; hardness ratios of spectra from events near the 
core can therefore be compared with hardness ratios at 
larger radii where pileup is known not to occur. This ef- 
fect is masked to some extent by the energy dependence 
of the Chandra PSF, which becomes broader as energy 
increases, so a more effective means of detecting pilcup 
is to study the distribution of events that are removed 
by VF filtering. 

In ACIS Faint (F) mode data, the distribution of 
charge within a 3 x 3-pixel 'island' is measured to de- 
termine if a valid X-ray event has been detected. VF 
filtering uses 5 x 5-pixel islands and discards events that 
have significant charge in the outer pixels. It is there- 
fore much more sensitive to event pilcup than F mode, 
and if an event is not removed by VF filtering it is al- 
most guaranteed to be a single unpiled event. We studied 
the distribution of events discarded by VF filtering as a 
function of radius for each observation phase and chose 
the excluded-core radii listed in Table Q] such that no 
more than 1.3% of the total flare-phase events outside 
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TABLE 1 

ACIS Spectral Data Extraction Parameters 



Gross Excluded Uuter btreak 

Time Exposure a Core Radius Boundary Width c 

Phase (2002 Sep 9 UT) (s) (pixels) (pixels) (pixels) 

Quiescent 00:01:20 - 11:51:40 4262(5 4 40 (radius) 4 

Flare 11:51:40 - 12:38:20 2800 8 202 (box height) 6 

Decay 12:38:20 - 16:51:45 15205 5 80 (radius) 5 

a The effective area corrections we apply during spectral analysis also account for detector 
deadtime (i.e., for exclusion of streak events). b The source core is excluded because of 
pilcup; core radii are listed for annular extractions. Streak extractions excluded pixels within 
20 pixels of the source to avoid contamination by annular events. c Streak data are only 
included for temporal analyses, with the listed box widths. Spectral extractions exclude a 
3-pixcl-wide box around the streak. 



that radius (0.7% for the quiescent phase and 0.8% for 
the decay) would be excluded by VF filtering. The net 
pileup fraction of the extracted events is estimated to be 
less than 0.5%. Close to the core but still within the ex- 
traction region, nearly all the 'bad VF' events are in fact 
valid and unpiled. To keep those events, and since the 
number of background events is negligible within such a 
small region, VF filtering was not applied out to a radius 
of twice the excluded core radius. 

The outer boundaries of the spectral extraction regions 
for each phase were chosen to keep the fraction of back- 
ground events small (< 4% for the Quiescent phase and 
less for the others) while including as many X-ray source 
events as possible. The readout streaks, which do not 
suffer from pileup because of their very short effective 
exposure times, were included in temporal analyses but 
excluded from spectral analyses (see Fig. [T] and Table [1]) 
because of concerns that the effective detector gain might 
be sufficiently different from that of non-streak events to 
distort the spectra. A g a in inc rease of ~7% has been 
reported by iHeinke et alj ([2003) in streak spectra from 
the front-illuminated ACIS-I3 chip at the Ir-M absorp- 
tion edge (a spectral feature near 2 keV arising from the 
Chandra mirrors), and an increase of ^2.5% was mea- 
sured in the back- illuminated ACIS-S3 chip 4 . It is not 
yet clear if this effect is energy dependent (private com- 
munication, R. J. Edgar), and our streak spectra had too 
few counts to measure a gain shift, other than to say that 
it must be no more than a few percent at the Ir-M edge 
and no more than 10% at energies down to ^500 cV. 

Response matrices (RMFs) and effective area files 
(ARFs) were created for each extraction region according 
to the psextract script, which uses the mkacisrmf tool. 
ARF correction functions were then created to account 
for exclusion of the source core as described below. 

4.2.1. Effective Area Corrections 

Although exclusion of the source core eliminates the 
spectrum-distorting effects of pilcup, which arc so severe 
in this case that they cannot be modeled, it introduces 
other modifications to the source spectrum because of the 
energy dependence of the telescope point spread function 
(PSF); low energy X-rays are more tightly focused than 
those with high energy. 

We modeled the PSF using the Chandra Ray Tracer 
(ChaRT; Carter et al. 2003) tool using 22 million rays 
and then projected them onto the ACIS detector us- 

4 See T. J. Gaetz (2004), at 

http : //exc .harvard. edu/cal/Hrma/psf /wing_analy sis .ps 



ing the MARX 5 simulator, producing 7.5 million de- 
tected events ranging from 0.1 to 8 keV with half of 
them below 3 keV. The MARX parameter DitherBlur 
was changed from the default value of 0.35" to 0.28" 
to properly model the effects of turning pixel random- 
ization off during reprocessing of the observation with 
acis_process_events. The data extraction regions 
listed in Table [T] were then applied to the simulated data, 
and extracted spectra were compared to the full-detector 
spectrum in order to derive effective area correction func- 
tions for each phase (see Fig. |4ji) . Near the peak of the 
detected source spectra around 1 keV, only 2 or 3% of 
the potential event detections are extracted. The average 
Quiescent count rate after correcting for the extraction 
efficiency is 2.0 cts s _1 . 

Fig-Hk also shows the extraction efficiency for the read- 
out streak during the Flare phase, using a 6-pixel- wide 
box and excluding a 20-pixel-radius around the source 
core. Readout-streak efficiencies are very similar for all 
three phases and nearly flat as a function of energy since 
the streak primarily comprises events from the core of 
the source, where energy dependent PSF effects are rel- 
atively minor. Some adjustments to the model streak 
results were required because MARX simulates full-chip 
(1024-row, 3.2-s frametime) ACIS operation while the 
observation used 206 rows with a 0.6-s frametime. The 
effective exposure time for the full-length streak is 7% 
larger in the observation than in the simulation 6 but less 
of the streak is extracted (162 of 206 rows versus 988 
of 1024, a relative difference of 23%), yielding a net ad- 
justment of ~14% to the extraction efficiency (observed 
lower than simulated). 

ChaRT simulations of the PSF are known to have small 
errors, and when excluding the core such errors may be 
relatively large compared to the remaining fraction of 
source flux. We therefore compared the simulation's pre- 
dictions of the ratio of annular and streak events versus 
the observed ratio for each observation phase. As seen in 
Fig. [4Jd, c, and d, uncertainties are dominated by observa- 
tional statistics and are especially large at high energies 
where there are few counts, but the simulations match 

5 |http : //space . mit' edu/CXC/MARX/ 

6 The streak is an artifact of the CCD readout process, with 
a net exposure time per frame equal to the number of CCD rows 
(206) multiplied by the time required to shift the image by one 
row during readout (40 /jsec), or 0.00824 sec. The readout-streak 
exposure efficiency when using a 0.6-sec frame-time is therefore 
0.00824/(0.6 + 0.00824) = 1.35%, and pilcup is completely neg- 
ligible. With the full array, 3.2 readout, the streak fraction is 
(1024 x (40 x 10 _6 ))/(3.2 + (1024 x (40 x 10~ 6 ))) = 1.26%. 
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Fig. 4. — Simulated and observed results on data extraction 
efficiency. Bins are 100 eV below 1 keV, and logarithmic (25 per 
decade) above 1 keV. (a): Only a small fraction of the total source 
flux is extracted using annular regions that exclude the source core, 
particularly at low energies where focusing is best. Extractions 
along the readout streak (excluding the core) have much less en- 
ergy dependence. (b,c,d): Ratio of counts in annular and streak 
extraction regions for Quiet, Flare, and Decay phases, for simu- 
lated and observed data, (e): Ratio of observed and simulated 
results from (b,c,d) for each phase. The observed ratio is larger 
than the simulated ratio at low energies, suggesting the simulated 
annular extraction efficiencies are too low. (/): Adjusted annular 
extraction efficiencies, after applying the smooth correction func- 
tions in (e) to results in (a). 



observation reasonably well except at low energies. These 
discrepancies are presumably due to underprediction of 
the enclosed count fraction in the outer core of the PSF 
at low energies where focusing is best. As noted in £|4.2I 
there may be small gain shifts in the streak spectra, but 
the effects on apparent streak extraction efficiency would 
be much smaller than the discrepancies seen. Further ev- 
idence that the simulation's effective area corrections are 
suspect below ^600 cV is provided by varying the energy 
range of spectral fits, as explained in 

Based on the data points in Fig.[4p-d, Fig.|4^ plots the 
ratio of simulated and observational results, along with 
smooth curves that approximately follow the results for 
each phase. Each curve shows how the ratio of simulated 
extraction efficiencies for annular and streak data regions 
would have to be adjusted to bring it into accord with 
the observed ratio. If one assumes that the simulated 
streak extraction efficiency is correct, as one would ex- 
pect since the PSF (or rather, line spread function) of 
the readout streak has so little energy dependence, then 
the smooth curve can be used as a correction factor for 
the annular-region extraction efficiency derived from the 
ChaRT/MARX simulation. Fig. [4f shows the simulation 
extraction efficiencies from Fig. 0Ja. multiplied by those 
empirical correction factors, yielding what we will refer 
to as "adjusted extraction efficiencies" . Uncertainties in 
the adjusted efficiencies are relatively large at higher en- 
ergies, but this is because of the small number of counts 
and the effect on spectral fits is modest. 

5. SPECTRAL ANALYSIS 

Spectral ana lysis was performed w ith the CIAO Sherpa 
fitting engine (|Freeman et al.ll2001h . First, the adjusted 
extraction efficiency curves (Fig. |4f) were fitted with 
smooth functions (the sum of five Gaussians) to be used 
in conjunction with the ARFs (as multiplicative factors) 
during spectral fitting. The coronal emission itself was 
modeled using xsvap ec models based on the plasma 
emission code APEC ([Smith et all l2001al ). Parameter 
estimation was performed using the modified \ 2 statistic 
(|Gehrelslll986t ). with some rebinning at higher energies 
to ensure at least 20 counts per bin. 

Fits extended down to 350 eV, below which the ACIS- 
S response becomes increasingly uncertain, with upper 
energy limits for each phase determined by the level of 
source emission versus background (see Table [2J • Differ- 
ent energy ranges were also tried, with variable effects 
on the fitting results. Fits of the Quiescent spectrum 
down to 200 eV led to overpredictions of ~20% in the 
flux below 500 cV, while fits using a lower limit of 500 
eV yielded poor constraints on the O abundance. Fits 
to the Flare and Decay spectra were less sensitive to the 
lower limit, and in all cases varying the upper limit of 
the fit range had relatively little effect. 

Because of the modest resolution and statistical qual- 
ity of our spectra, we grouped elements with similar 
first ionization potential (FIP; in parentheses) as fol- 
lows (see §5.41 for the reasoning behind FIP grouping): 
C, N, O, and S (11.3, 14.5, 13.6, and 10.4 eV); Na, 
Mg, Al, and Si (5.1, 7.6, 6.0, 8.2 eV); Ca, Fe, and Ni 
(6.1, 7.4, 7.6 eV); Ne and Ar (21.6, 15.8 eV). Element 
abundances were dete rmined relative to the solar photo- 
spheric abundances of I Anders fe Grevessd (|1989f ). which 
is the Sherpa default. More recent photospheric abun- 
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TABLE 2 
Spectral Fitting Results 



Quiescent ~ Flare" Decay 3 " 



Parameter 


1-T 


2-1 


1-T 


'1-T 


'1-T 


Fit range (keV) 


0.35- 


3.0 


0.35- 


8.0 


0.35-8.0 


Counts in fit 


1802 


4829 


4891 


Est. background 


74 


65 


172 


Deg. of freedom 


63 


61 


164 


162 


125 


Reduced x 2 


0.61 


0.53 


0.99 


0.78 


0.71 


fcTi (keV) 


46+ 03 


0.98^ 


2 95+ - 27 
z - 3J -0 .20 


o 7 i +4.30 
J - ' -0.26 


1 91+012 
^ -0.12 


kT 2 (keV) 




0.29±°;j» 




q 41+1-17 


n 37+0.04 


EMi b (10 50 cm" 3 ) 




4 1+5.8 


o 4 i+36 


971 

z,1 -20 


1 o 7+2.8 
lz -'-2.4 


EM 2 b (10 50 cm" 3 ) 








26l 306 


97 O+3.0 
z ' ,o -3.0 


Lum. c (10 27 erg s^ 1 ) 


9.35 


8.66 


449 


464 


51.1 


OCNS abundance d 


12+ 005 

u ' lz -0.04 


13+ 007 


21+ ' 47 
u - zl -0.21 


n ,0+0.26 
U-3»_ 20 


4O+0.17 
u -^ z -0.13 


NeAr abundancc d 


o 2q+ 07 


34+ - 24 


1 39+0-90 


z,lo -0.59 


1 fi^+O- 48 


MgNaAISi abund. d 


10+ 006 
u - lu -0.06 


15+ 020 
U ' i0 -0.13 


u ' uu -0.00 


44+ - 54 


"Sfi+0- 32 
U.00_ 25 


FeCaNi abund. d 


04+ 01 


1Q+ 010 


27+ ' 14 

U - Z '-0.13 


41+O.I8 
-0.38 


40+° 13 



Note. — Quoted uncertainties are formal 68% confidence intervals for the fits and 
do not include systematic uncertainties, such as those for the effective area of a data 
extraction region. 

a Underlying quiescent emission (1-T model) is included in the flare and decay 
fits. b Emission measure is defined as J n e riH dV and is equal to the Sherpa fit 
normalization (adjusted for pileup) xlO 14 X 4nD 2 where D is the source distance (2.97 
pc) in cm. c Luminosities (averages over each phase) are for the 0.25-11 keV band. Flare 
and tail luminosities are in addition to the underlying quiescent emission. d Dominant 
clement within each grouping is listed first. Abundances are relative to solar photosphcric 
values listed by Anders Sz Grevesse 1989. 

dance tabulations are available (e.g.. iGrevesse fc Sauvall 
199a iAsplund~G revesse. fe Sauval 1 120051) . but the effect 
on our fit results (in terms of abundances relative to H) of 
using different assumptions is negligible, as confirmed by 
trials using the Asplund et al. (2005) abundances, and as 
expected since each element grouping typically has one 
dominant species. O emission dominates that from N and 
C because the ACIS effective area decreases rapidly to- 
ward lower energies, and there is virtually no S emission 
in the quiescent spectrum. Ne emission similarly domi- 
nates that from Ar. Abundances for Mg and Si, which 
have very similar FIP, arc roughly an order of magnitude 
larger than those for Na and Al, and Fe likewise dom- 
inates Ca and Ni. Abundance linkages were studied in 
more detail for the flare spectrum, as described in ij5.2l 
Our analysis assumes that quiescent-phase emission 
is always present at a constant level, with added com- 
ponents for the flare and decay phases representing lo- 
calized emission. Quiescent fit parameters are therefore 
frozen during flare and decay fitting, with the normaliza- 
tion adjusted for the differing exposure time and extrac- 
tion region efficiency of each phase. The quiescent emis- 
sion is completely overshadowed by higher-temperature 
components during the flare and decay, however, so its 
inclusion has little effect on those fits. The background is 
scaled and subtracted and contributes ~4% of the counts 
in the fits' energy ranges (less for the flare). Abun- 
dances are free to vary for each phase, subject to the 
FIP grouping described above. Interstellar absorption to 
this nearby source is completely negligible at the energies 
of relevance here (Nh < 1 x 10 18 cm 2 ; Wood et al. 2005). 
Fit results are summarized in Table [5] and described in 
detail below, and the spectra are shown in Figure [51 

5.1. Quiescent Spectrum 

As can be seen in Table [5J results for one- and two- 
temperature fits to the quiescent spectrum are quite sim- 



ilar, apart from a marginally significant difference in Fe 
abundances. Using two temperatures reduces the % 2 
slightly, but the F-tcst significance of the second com- 
ponent is 0.29, much larger than the typical threshold of 
0.05 for using a more complex model. We therefore refer 
to the 1-T fit in all subsequent discussions unless other- 
wise noted. The fact that the reduced \ 2 is significantly 
less than 1 (0.61) is likely a reflection of the modest num- 
ber of counts. The similarity of abundances from the 1-T 
and 1-T fits, however, suggests that the derived values 
are realistic. 

The total X-ray luminosity is ~ 9 x 10 27 erg s" 1 , some- 
what higher than the value o f 6 x 10 27 erg s" 1 est imated 
from ROSAT observations (jHunsch et al J 1 1999( 1. The 
difference is probably due to a combination of real source 
variability and uncertainties in spectral modeling; uncer- 
tainties in spectral extraction efficiencies (Fig. Hh) are 
largest the low energies, leading to relatively larger un- 
certainties in luminosity during quiescence than in the 
flare and decay phases. 

The best-fit temperature is T = 5.0 x 10 6 K (kT = 0.43 
keV) and the 68% confidence intervals for the abundances 
arc 0.29t°;° 7 for Nc (and Ar), 0.12±8;{g for O (and C, N, 
S), 0.10±g;^ for Mg (and Na, Al, Si), and 0.04tg$ for 
Fe (and Ca, Ni; the 1-T fit had higher Fe abundance but 
with large uncertainties). Abundance-FIP correlations 
are discussed in detail in §5.41 but we note that in qui- 
escence, the derived element abundances correlate (with 
modest significance) with FIP. 

To study relative abundances in more detail and aid 
comparisons of quiescent and flare abundances, we stud- 
ied x 2 as a function of the Ne/O abundance ratio by 
freezing Ne/O and refitting the data (see Figure [6]). The 
result is a best-fit ratio of 2.4 with a la confidence range 
(95.5%, Ax 2 = 4.0) for between 1.3 and 4.1, somewhat 
less than the value found from the flare spectrum (see 
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Fig. 5. — Fits to spectral data, using the Quiescent "1-T," Flare 
"2-T," and Decay models listed in Table [2] Data are background- 
subtracted. 



j)5.2[) . A study of the Fe/O ratio shows a similar en- 
hancement during the flare, but the significance of this 
result vanishes if the quiescent Fe abundance from the 
2-T fit is used instead of that from the 1-T fit. Uncer- 
tainties on the flare Mg abundance are too large to draw 
any conclusions regarding a difference from the quiescent 
abundance. 

5.2. Flare Spectrum 

One- and two-temperature models were also employed 
to fit the flare spectrum with, as noted above, an un- 
derlying fixed quiescent component. While the 1-T fit 
gives a formally adequate fit with reduced x 2 = 0.99, 
the 2-T model provides a visibly much better fit at low 
and high energies. We therefore used the 2-T model, al- 
though both fits yield consistent results with respect to 
element abundances. Note that the upper limits on the 
two components' temperatures are not well constrained, 
leading to a large uncertainty on the cooler component's 
emission measure. 

The power of this flare is remarkable. At its peak (ap- 
proximately 3.8 times the average power during the flare 
phase we defined) , the energy flux that would be observed 
1 AU from the star reached 0.34 W m~ 2 in the GOES 
solar flare energy band (1-8 A, or 1.55-12.4 keV), cor- 
responding to an X3400 solar flare. Over the full X-ray 
band (0.25-11 keV), the peak luminosity is ~ 1.8 x 10 30 
erg s" 1 (13% of the L^i listed by Fleming, Schmitt, & 
Giampapa 1995), and the total radiated energy including 
the extrapolated decay phase is ~ 2.3 x 10 33 erg. The 
only flares from isolated M stars significantly more ener- 
getic than this were observed on EV Lac by the Adv anced 
Satellite for Cosmology (A SCA: iFavata et al1l2000D and 
EQ1839.6+8002 by Ginga (jPan et al.lll997f) . Both those 
events were roughly ten times more powerful than the 
Ross 154 flare. 

In addition to its power and high temperature — note 
the Fe XXV emission feature at 6.7 keV in Figure [5] — an 
interesting aspect of the flare is the large enhancement 
of Ne relative to its quiescent value. While all the abun- 
dances rise during the flare (uncertainties for Mg are too 
large to draw any conclusions), the increase for Ne is 
the most significant. Because abundance ratios are often 
more reliable than absolute values, we again applied the 
analysis procedure described in H5.1l to study the signif- 
icance of the Ne/O ratio (see Figure [5]). We conclude 
that there is roughly a 90% statistical likelihood that Ne 
is relatively more abundant during the flare than during 
quiescence, with a best-fit Ne/O ratio of 5.7, compared 
to the quiescent ratio of 2.4. The 90% statistical like- 
lihood, however, docs not take into account systematic 
uncertainties from our extraction efficiency modeling. 

The flare was hot enough to generate substantial emis- 
sion from H-like and Hc-like Ar, unlike during the qui- 
escent phase, so we unlinked the Ar and Ne abundances 
in an alternative fit. The best-fit results yielded very 
large uncertainties for Ar with negligible effect on the 
Ne abundance or any other parameters. Fits that un- 
linked S (from O, C, and N) and Si (from Mg, Na, and 
Al) provided similarly uninformativc results. 

As noted above, the interstellar column density to Ross 
154 is so low (Njj < 10 18 cm" 2 ) that X-ray absorption is 
negligible. There has been one report, however, of a large 
increase in absorbing column density (to ~ 2.7 x 10 21 
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Fig. 6. — Confidence intervals for the Ne/O abundance ratio 
(normalized to solar photospheric abundances) in fits to quiescent 
(solid curves) and flare (dashes) spectra. The 1-T model was used 
for the quiescent fits and the 2-T model for the flare (see Table 2). 

cm~ 2 ) during a large flare observed on Algol (B8 V + 
K2 IV) by BeppoSAX, presumably the result of a coronal 
mass ejection fr om the K2 secondary in association with 
the flare onset (jFavata fc Schmittl Il999f ). We therefore 
added an absorption term to our flare model, but the fit 
drives Nh to zero. Freezing Nh has no significant effect 
on the fit for values up to a few times 10 20 cm -2 , but fits 
with Nh much beyond 10 21 cm -2 are noticeably poor. 

5.3. Decay Spectrum 

One-temperature models give unacceptable fits to the 
decay spectrum and Table [5] therefore lists results only 
from 2-T models. Physically, one would expect temper- 
atures and abundances to be intermediate between the 
quiescent and flare values, and probably closer to the lat- 
ter since we observe the portion of the long-lived decay 
immediately following the flare. The fit therefore gives 
sensible results, but the model and fit uncertainties are 
too large to draw any firm conclusions regarding clement 
depletions or enhancements. 

Before discussing implications of the fit results, it is 
worth recalling that the real uncertainty in any phys- 
ical parameter is always larger than the formal fitting 
error. Although the fit results in Table [2] are fairly in- 
sensitive to the choice of model, model assumptions in- 
evitably affect the values of derived parameters, as do 
many other factors such as atomic data uncertainties, 
limited energy resolution, and calibration errors. The 
latter can be especially important with regard to the 
interdependence of derived emissio n measures and ele- 
ment abundances , as illustrated bv iRobrade fc Schmitt) 
(|2005f ) in fits to XMM-Newton data from several in- 
struments; much of the (modest) disagreement be- 
tween results from fits to the MOS and other instru- 
ments can likely be attributed to recently discovered 
spatially dependent response function variations in the 
MOS detectors caused by cumulative rad iation exposure 
(|Stuhlinger et al.ll2006t iRead et al.ll2006f ). Although cal- 
ibration errors in ACIS should not be a significant is- 
sue here, exclusion of the piled-up source core and at- 
tendant energy-dependent modifications to the effective 
area have introduced difficult-to-estimate uncertainties 
into our fits. Uncertainties in atomic data used in the 
fits are relatively small and, as described above, we have 



taken some care to study the sensitivity of our model 
assumptions and parameters, but the uncertainties on 
the fit values listed in Table [2] should be viewed as lower 
limits. 

5.4. Abundances Discussion 

Coronal abundance patterns seen in stars with solar- 
like to very high activity levels currently appear to de- 
viate from expected photospheric values according to 
element first ionization potentials. In the solar case, 
low-FIP elements (< 10 eV, e.g., Si, Fe, and Mg) ap- 
pear enhanced by typical factors of 2-4 compared with 
high-FIP elements (> 10 cV, e.g, O, Ne, Ar) which 
have roughly photospheric values (e.g.. iFeldmanl U992L 
and references therein). In stars, there appears to be a 
steady transition from a solar-like FIP effect to "inverse- 
FIP effect" as activity level rises (e . g..lDrake et al.lll995l: 
White et all H99l iTeiTeschi et all 12001 "^ 



20031: IGudel eTaif 2002b) . 
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Several abundance studies of active M dwarfs suggest 
they s hare a pattern simila r to th at of hig her mass active 
stars. IRobrade fc Schmittl (|2005l . see also iRaassen et all 
120031 ) analyzed XMM-Newton observations of four M 
dwarfs with spectral types similar to that of Ross 154, 
namely EQ Peg (GJ 896AB; a wide binary with M3.5Vc 
and M4.5Ve components), AT Mic (GJ 799 AB; a tidally- 
interacting M4.5Ve+M4V binary), EV Lac (GJ 873; 
M3.5Ve), and AD Leo (GJ 388; M3.5e). They con- 
clude that all these stars exhibit a similar abundance pat- 
tern: a remarkably flat abundance-versus-FIP relation- 
ship (roughly half-solar photospheric) from Si (8.15 eV) 
to N (14.53 eV), with only Ne (21.56 eV) showing a 
conspicuous enhancement. This pattern is also rem- 
iniscent of t hat fo und for Proxima Cen (GJ 551) by 
IGudel et alj (|2004l ) and for the dMe active spectro scopic 
binary YY Gem (Castor C) bv IGudel et "all (|2001h from 
XMM-Newton observations. 

Adopting the 2-T abundance results in Table [21 we 
find Ross 154 is broadly consistent in relative abundances 
with those from the earlier studies of M dwarfs. We find 
no significant deviation from a solar mixture in terms of 
the relative abundances of the OCNS, MgNaAISi, and 
FeCaNi groups, but a significantly larger relative abun- 
dance of the NcAr group in quiescent, flare and decay 
spectra. There is some indication that the coronal metal- 
licity (i.e. all metals scaled together) differs between qui- 
escent and flare/decay phases: both FeCaNi and OCNS 
groups appear more consis t ent w ith a value of 0.4 times 
that of lAnders fc Grevessd (jl989t ) during the decay, com- 
pared with ~ 0.15 during quiescence. This is confirmed 
by 2-T model fits in which all elements were tied to their 
relative solar values (except for the Ne/O ratio which was 
fixed at 2.4) and only the overall metallicity was allowed 
to vary. The metallicity rises from the quiescent value 
of ~ 0.15 to ~ 0.55 in the flare, with a best-fit value of 
~ 0.45 during the decay (but statistically consistent with 
the flare value). The flare and decay metal abundances 
are essentiall y the same as the photospheric estimate of 
lEggenl (|l996f ). 

The Ross 154 abundance patter n and Ne/O abun- 
dance ratio is similar to those of the IRobrade fc Schmittl 
(|2005| ) sample, although during the flare there is a 
suggestion at the 1.5<7-2cr level of a Ne/O increase 
(see and Tabic |2|). The Ne/O ratio is inter- 
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csting in the context of understanding solar structure 
and its distribution in the coronae sta r s of different 
activ i ty level (e.g., Basu fe Antial 120041; iBahcall et al.1 



\2Ml IDrake fc Testall2005l) . IDrake fe Testal (120051) found 



a remarkably constant Ne/O abundance ratio that 
is ~2.4 times higher than the currently favored so- 



lar ratio of 0.15 by number (e.g., 


Anders & Grevesse 


19891: lAsplund, Grevesse, & Sauval 1 


20051: Younal 2005: 


Schmelz et all 120051: iLandi, Fcldman, & Doschekl 12007) 



in a sample of mostly active stars over a wide range of 
spectral type (and as seen in Ross 154). They argue that 
a FlP-based fractionation mechanism is unlikely to re- 
sult in such consistent Ne/O ratios when the Fe/O ratio 
in th e same stars varies over an order of magnitude (see, 
e.g., [Gudcl 20041). and suggest instead that the higher 
ratio — including tha t found for the M dwarfs analyzed by 
iRobrade fc Schmittl (|2005l ) — represents underlying pho- 
tosphcric values. In this scenario, the M dwarfs appear 
to exhibit coronal abundances that reflect their relative 
photosphcric values but that differ in terms of a global 
metal depiction factor. 

The quiescent cor onal metallicity of Ross 154 is some- 
what lower than the IRobrade fc Schmittl (|2005[ ) sample, 
and the difference in quiescent activity level between 
Ross 154 and the other dMe's may be important here. 
The X-ray luminosities of stars in the latter group range 
from a few 10 28 erg s _1 up to several 10 29 erg s , on av- 
erage roughly a factor of 10 higher than that of Ross 154. 
Ross 154 might then be probing a lower activity regime 
that has a characteristically larger coronal metal deple- 
tion. At face value, it would appear from Ross 154 that 
the degree of depletion increases with decreasing activity. 
However, Proxima Cen is a star similar to Ross 154 in 
size but with somewhat lower activity (Lx ~ 2 x 10 27 erg 
s _1 ) and an essentially photosp heric (solar- like) coronal 
composition (|Giidel et alJ 12004 ). Similarly, YY Gem is 
a very active system but appears to hav e sligh tly lower 
abundances than the IRobrade fc Schmittl (|2005| ) sample. 
Both stars thus suggest the opposite trend: increasing 
coronal metal depletion with increasing activity, as is ob- 
served for higher mass dwarfs. 

What, then, is the underlying explanation for 
the different M dwarf coronal abundances? The 
IRobrade fc Schmittl 1)20051 ) sample includes active single 
stars that must be relatively young and composition- 
ally representative of the local cosmos. Since the ob- 
served coronal metal abundances a factor of two lower 
than the local cosmic (essentially solar) values, it seems 
likely that the more active M dwarf coronae are de- 
pleted in metals relative to photospheric values by fac- 
tors of ~ 2. This view is supported by the increase in 
coronal metals seen in the Ross 154 flare: such a large 
flare probably introduced fresh material from the chro- 
mosphere or photosphere into the corona via chromo- 
spheric evaporation, as conjectured to explain similar 
abundance changes seen durin g large flares on other M 
dwar fs (e.g. iFavata et al.ll2000l) and interacting binaries 
(e.g.. lFavata fc Schmittl fl999l : iGudel et al.lll999h . Other 
differences from the expected coronal metallicity-activity 
trend are likely to be present simply due to scatter in 
photospheric metallicity of the young, active stella r pop- 
ulation (e.g.. lBensbv. Feltzing. fc Lundstr5ml l2003). The 
slightly s ub-solar photo spheric metallicity of Ross 154 
found by lEggenl ()1996l ) supports this conjecture: its 



corona is metal poor as a result of both low photospheric 
metallicity and some degree of coronal metal depletion. 

Such an abundance pattern is difficult to confirm in 
M dwarfs. There are few low- activity/low- luminosity 
dMe stars nearby enough for detailed study and thus 
nearly all well-studied dMe stars are highly active. Nev- 
ertheless, these intriguing results suggest that more de- 
tailed observations of Ross 154 and other more inactive 
stars would be highly worthwhile. 

6. QUIESCENT MICROFLARING 

Although the origin of coronal heating remains one of 
the great unsolved problems of modern astrophysics, it 
has been well established that solar flares are distributed 
as a power law in energy with the form ^ oc E~ a 
dLin et all fl98i iHudsonl Il99lh . The index a is ap- 
proximately 1.8 for high-energy flares, possibly increas- 
ing to ~2.6 for low-ener gy flares (|Krucker fc Benj|1998l : 
IWinebarger et al.1 I2002D. although this latter m easure- 
ment is disputed ( Aschwanden fc Parnelll I2002D . The 
precise value of a is of fundamental importance because 
if a > 2 then it is theoretically possible to attribute the 
entire coronal heating budget to energy deposited into 
the corona during flares. 

Several authors (iButler et al.l 119861: 



Ambruster. Sciortino fc Golubl 11987 : Robinson et al.l 



19951: lAudard et all 120001 : (Kashvap et all [2003 



Giidel et all 120031 . I2004D have therefore suggested 



that the apparently quiescent emission in other stellar 
coronae may be due to the continuous eruption of small 
flares, commonly referred to as "microflaring," although 
the flares involved extend well into the range of solar 
X-class events. The luminosity of these microflares is 



roughly 10 -10 ergs s ; for comparison, a solar XI 
flare corresponds to a few times 10 26 ergs s _1 over the 
full X-ray band. Currently, the detectability limit for 
flares on stars other than the Sun is a l uminosity of 
roughl y 2 x 10 26 ergs s" 1 (for Proxima Cen: IGudel et al.1 
I2002al) . 

For large flares observed with EUVE, Audard et al. 
(2000) found that the majority of cool stars they ana- 
lyzed had a > 2. Using a more sophisticated method 
that directly analyzes photon arrival times and thus ex- 
tends quantitative modeling to weaker flares, Kashyap 
et al. (2002) found that low-mass active stars such as 
FKAqr, V1054Oph, and AD Leo all have a > 2, and 
that a large fraction of the observed emission (generally 
> 50%) can be attributed to the influence of numer- 
ous small overlapping flares. Their method analyzes the 
distribution of arrival time distributions that would be 
realized for different flare distributions, with the flares 
themselves modeled as randomly occuring in time with 
instantaneous rises and exponential decays. The intensi- 
ties of individual flares are sampled from the aforemen- 
tioned power law. 

Because the model is stochastic in nature, Monte Carlo 
methods are used to obtain best-fit values and errors of 
the flare intensity and the power-law slope index a. For 
a given X-ray microflaring luminosity, smaller values of 
a yield light curves dominated by occasional large flares, 
whereas light curves characterized by larger values of a 
are dominated by numerous smaller flares. This differ- 
ence in the general characteristics of light curves can be 
exploited to determine how the flares arc distributed on 
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the star. Naturally, better constraints are obtained on 
a when photon arrival times are known to higher accu- 
racy, especially for smaller values of a since such models 
contain flares with larger peak count rates and therefore 
less time between events. 

6.1. ACIS Observation Results 

The analyses described above model flares with single- 
exponential decays. We use the same flare model with 
the decay time constant set to 1000 s, a value character- 
istic of flares at typical coronal densities; the actual value 
of the decay time has little effect on our results for the 
relatively small flares analyzed here. Because the large 
flare in the ACIS observation clearly decays with double- 
exponential behavior (see ^Z]), we restrict our analysis to 
the quiescent phase. 

The resultant two-dimensional joint posterior proba- 
bility density p(a,rp) of the parameters rp, the event 
rate due to flares, and a, the flare distribution index, is 
shown in Fig. [7] The 1-dimcnsional probability density 
on a can be obtained by integrating over rp, 7 and we 
find that a lies between 1.96 and 2.86 (90% confidence 
level) with a mean of 2.45. Note that p(a, rp) is not well 
localized, and the flare distribution is only poorly charac- 
terized by crude summaries such as the means and vari- 
ances of tf and a. The figure shows that smaller values 
of a (corresponding to relatively fewer but larger flares) 
are consistent with a smaller rp, i.e., the contribution of 
flares to the total emission is smaller. In contrast, for 
larger a, when the light curve is dominated by more fre- 
quent but less energetic flares, the contribution of flares 
to the total count rate becomes more significant. 

6.2. HRC Observation Results 

The complex probability distribution from the ACIS 
analysis is revealed more clearly by the HRC data, which 
have more counts (20500 X-ray events plus ^10000 back- 
ground events in 46776 s for the HRC observation, ex- 
cluding flares; ~ 4400 X-ray events in 42620 s for the qui- 
escent ACIS data) and better temporal resolution (^3.5 
ms versus 0.6 s). The relationship between temporal res- 
olution and the accuracy of a is complex and has not 
been investigated here because of the extreme compu- 
tational demands of such an analysis, but some of the 
differences between the ACIS and HRC results discussed 
below are likely due to differences in the accuracy of event 
timing. In particular, the HRC's superior timing permits 
better characterization of p(a, rp) at smaller values of a. 

As described in Sj3l the HRC-I data were collected far 
off-axis, so that the source image was spread out over 
a large area. Detector background was steady at 2.55 
cts s _1 over the full field (1.86 cts s" 1 after standard 
filtering) and the background-subtracted quiescent X-ray 
rate was 0.3-0.5 cts s _1 (see Fig. [5]). A large source 
extraction region (an ellipse with major and minor axes 
of 400 and 600 pixels; see Fig. [3]) was required to collect 
the X-ray events, but because the background is steady, 
the microflaring analysis is not sensitive to background 

7 The joint posterior probability density function (pdf) is techni- 
cally written as p(a,rp\D), where the notation indicates the prob- 
ability of the parameter values given the data, D. For the sake of 
notational simplicity, we drop the explicit conditional in all refer- 
ences to the posterior pdf. Thus, the posterior pdf of a marginal- 
ized over rp is referred to as p(a) = p(a\D) = J p(a,rp\D) drp. 



contamination. Indeed, the events are modeled with a 
constant non-flaring base event rate that includes the 
background rate. 

A large flare occurred near the very end of the obser- 
vation, exceeding the telemetry limit of ~3.5 cts s _1 and 
leading to a loss of events despite some buffering capabil- 
ity. Because of the limitations of the Next-In-Line mode, 
the deadtime can not be determined directly, but the 
flare intensity was at least a factor of 4 higher than the 
quiescent rate and probably close to a factor of 10 (see 
below). Some lesser flaring at 3 or more times the qui- 
escent rate was observed earlier in the observation, and 
a coincidental dip in the telemetered background rate 
indicates that the total event rate exceeded telemetry 
capacity during that time. Telemetry buffering appears 
to have preserved all the events during a few prior peri- 
ods (between ~ t = 12 and 14 ks) when the rate briefly 
exceeded 3.5 cts s _1 . 

Note that the upper curve in Fig. [8] plots the event 
rate for unfiltcred full-field data, which is most relevant 
for comparisons with the telemetry limit. The source and 
background light curves are for filtered data. Standard 
filtering reduces the HRC-I background by 25-30% with 
a loss of ~2% of valid X-ray events. Assuming that the 
true background rate is constant, one can derive the de- 
tector livetime during the flares by dividing the measured 
background rate by the true rate (derived from non-flare 
periods). The red traces in Fig. [5] plot the deadtimc- 
corrected source rate during the two flares, with a statis- 
tical accuracy of ~10% for the 100-s bins. Because many 
events during the large flare were lost, and those that 
remain may have incorrect times, we exclude the flare 
from our microflaring analysis. Lost events during the 
moderate flare near the middle of the observation may 
have affected the results so we repeated the analysis after 
excluding 1.4 ks around t ~ 17 ks. 

In both cases, we see in Fig. [7] that the HRC-derived 
joint probability distribution p(a, rp) is grossly similar to 
that determined from the ACIS data, but there are some 
important differences. First, at the 90% confidence level, 
we can rule out a < 2, since we find that 1.99 < a < 2.84 
(when excluding the large flare) and 2.08 < a < 2.87 
(excluding both flares). More importantly, given that 
p{a) provides a rather poor summary of the complex 
2-D p(a,rp) distribution, the broad correlation found 
between rp and a with ACIS data is now resolved into 
two major components in the probability distribution: 
one with a m 1.8, similar to the solar case, corresponding 
to a low value of rp < 0.1 ct s~ x , and another where the 
observed source counts are almost entirely dominated by 
flaring plasma, with rp m 0.4 ct s _1 , and a m 2.5. 

This bimodal distribution is an indication that the 
model adopted for fitting, that of flare intensities dis- 
tributed as a power-law with a single index, is too sim- 
plistic. The true flare distribution may be a broken 
power-law or another more complex form that requires 
consideration of local plasma conditions such as density, 
composition, and emissivity. For example, the correla- 
tion between energy deposition and X-ray emission may 
break down for small flares. These data thus provide the 
first concrete hint that different types of processes occur 
on late type stars than occur on the Sun, where the exis- 
tence of a single power-law spanning at least 6 orders of 
magnitude in flare energies (from 10 26 to 10 32 ergs) has 
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Fig. 7. — Joint probability distributions p(a,rp) of the power- 
law flare-activity index a and the average contribution of flares 
to the observed event rate, rp, for quiescent emission during the 
ACIS-S and HRC-I observations. The horizontal (blue) line near 
the top of each plot marks the total quiescent event rate with its 
lcr uncertainty range denoted by the dashed lines. The complex 
contours below that enclose 68%, 90%, and 95% probabilities for 
p(a, rp), with the large cross locating the mean values of a and rp 
and the arms signifying the lcr confidence ranges. Although it can 
be seen that the probability distributions are more complex than 
simple ellipses, on average roughly half the total quiescent emission 
can be ascribed to microflaring based on the ratio of rp and rtot- 
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Fig. 8. — Light curves of the HRC observation, with 100-s bin- 
ning. Zero time corresponds to the beginning of the observation on 
2007 May 28 at 4:19:18 UT. Thin vertical lines mark flare intervals, 
and the red traces within those intervals denote event rates after 
approximate corrections for deadtime due to telemetry saturation. 
The lowest curve is from the background region in Fig. E] and is 
scaled for the area of the source region. The microflaring analysis 
excludes the large flare at the end (after time=46842 s). Exclusion 
of the small flare between times 16642 and 18042 s has little effect 
on the analysis. 

been well established ([Aschwanden et aLll2000f ) . 

7. TEMPORAL ANALYSIS OF THE ACIS FLARE 

Having completed discussion of the HRC observation, 
we return to the ACIS data. As seen in Fig. ® the X-ray 
count rate during the flare increased by over a factor of 
100 from its quiescent level. The increase in luminosity 
is a factor of ~ 200, and as noted in £15. 2( observations 
of such large flares are quite rare. Observations of dMe 
stars with simultaneous X-ray and optical coverage are 
also uncommon but include: BY Dra with EXOSAT (de 
Jager et al. 1986); UV Ceti with EXOSAT and ROSAT 
(Butler et al. 1986; Schmitt, Haisch, & Barwig 1993); 
EQ Peg with EXOSAT and ROSAT (Butler et al. 1986; 
Katsova, Livshits, & Schmitt 2002); Proxima Cen with 
XMM-Newton (Gudel et al. 2002a); and EV Lac with 
Chandra (Osten ct al. 2005). 

In the two-ribbon solar flare model, magnetic recon- 
nection energizes electrons in the corona that are then 
accelerated into the chromosphere where they heat and 
ionize plasma that explosively expands into the corona in 
the process of "chromosphcric evaporation." The multi- 
MK plasma then cools as it emits soft X-rays. At the 
beginning of the flare, synchrotron radio emission is pro- 
duced as electrons spiral along the magnetic field lines 
and nonthermal bremsstrahlung is emitted as the elec- 
trons crash into the chromosphere at the magnetic field 
footprints. Optical emission is closely correlated with the 
nonth ermal hard X-ray (E > 10 keV) emission in solar 
flares (jHudson et al.lfl992T l and is often used as a proxy 
for the latter. Nonthermal flare emission is orders of 
magnitude weaker than the soft X-ray thermal emission 
and has been observed tentatively only once in another 
star, during a sup erflare in the active binary II Pegasi 
(|Osten et al.ll2007l ). 

Soft and hard X- ray light curv es are often related by 
the Neupert effect (|Neupertlll968l ). Nonthermal hard X- 
ray (and optical) emission traces the rate of electron en- 
ergy deposition in the chromosphere, while the soft X-ray 
emission is approximately proportional to the cumulative 
thermal energy transferred to the flare plasma by the 
electrons, at least during the beginning of the flare be- 



12 



Wargelin et al. 



Data (0.25-8 keV) 

Sum of exponentials 

Fast decay (t=400 s) 

Slow decay (t= 13000 s) 




Fig. 9. — X-ray light curve of the flare and its decay (top) and 
detail of the flare in different energy ranges (bottom); a slightly 
less restrictive spatial filter was used to include more counts for 
the bottom panel. Bin size is 100 s (20.5 s for optical in bottom 
panel). The pre-flare quiescent count-rate average has been sub- 
tracted from the optical data. The three data points near t = 43 
ks with wide bars in the optical curve represent average values over 
the indicated time intervals; gaps are dither-induced dropouts (see 
i|4.1H . The bottom-most curve traces the time derivative of the full- 
range (0.25—8 keV) X-ray light curve during the main flare. Beyond 
t = 44.2 ks (the flare peak) the derivative is negative. Prior to the 
main flare the statistics are poor, but the derivative light curve 
seems to peak roughly in concert with the initial optical flare at 
t m 42.7 ks, as one would expect for the Neupert effect. The 1-19 
bar across the top of the bottom panel refers to the 150-s bins used 
in Figs. rUTl and ITT1 



fore significant cooling has occurred. The time derivative 
of the soft X-ray light curve is therefore proportional to 
the hard X-ray light curve. Originally observed in solar 
flares, the Neupert effect (using optical or radio emission 
in lieu of hard X-rays) h as also been observe d in a few 
stars, including AD Leo (jHawlev et al.lfl995i ). Proxima 
Cen (jGiidel et al.ll2002af ). and probably II Peg (directly 
in hard X-rays; Osten et al. 2007). 

The Ross 154 flare has two subphases: an initial count- 
rate rise of a factor of ten over a period of ~500 s fol- 
lowed by a ~500-s plateau, and then the main flare with 
another ten-fold rise. Both the initial and main flares 
are accompanied by optical flares (see Fig. [9] and note 
that the pre-flare quiescent average of 28446 ± 264 cts 
bin -1 has been subtracted from the optical data). In the 



initial flare, the optical signal rises abruptly to its peak 
(within ~50 s) while the X-ray signal continues to rise 
after the optical flare has started to decline. Although 
the statistics are limited, this behavior is consistent with 
the Neupert effect described above, in which the optical 
signal is proportional to the time derivative of the X-ray 
signal until the X-ray peak is reached and plasma cooling 
takes over. 

In the main flare, the optical and X-ray light curves 
roughly mirror each other. This is particularly true for 
the harder X-rays (above 2 keV), which peak ~100 s be- 
fore the softer X-rays, in concert with the optical light 
curve. The X-ray and optical light curves also begin 
their rise at about the same time, but the optical signal 
increases more quickly at first and then slows its rise to 
the peak. In this case, there is no simple Neupert ef- 
fect, as the optical emission continues to rise even after 
the X-ray derivative has started falling. One possibility 
is that magnetic reconnection continues after the initial 
impulse, or perhaps conductive heating of the chromo- 
sphere by the flare maintains optical emission while the 
flare decays. A significant minority of large solar flares do 
not show a Neupert effect, and a similar mix of flare be- 
haviors is seen on other stars (e.g., Proxima Cen; Giidcl 
et al. 2004). 

7.1. Size and Density Estimates 

Although detailed dynamical and structural modeling 
of the flare on Ross 154 is beyond the scope of the present 
paper, simple geometric arguments can be used to derive 
rough estimates of the size of the flaring volume and the 
plasma density within it. We focus on the main flare 
event, which occurs on three time scales. First is the rise 
of the flare (r r i se ~ 500 s) followed by its decline, which 
can be modeled as the sum of a fast (r/ QSt « 400 s e- 
folding time) and slow (t s i ow « 13000 s) exponential de- 
cay. This double-exponential decay is often seen in stellar 
flares (e.g., Osten & Brown 1999; Giidel et al. 2004; Reale 
et al. 2004) and is also similar to the behavior seen in 
large solar flares, in which a steady active region under- 
goes a strong magnetic reconnection event resulting in an 
intense flare, and is followed by an arcade of reconnected 
loops that slowly decay (e.g., th e Bastille Day flare, see 
lAschwanden fc Alexand er 2001). Such a scenario is also 
supported by the temporal analysis of §E\ that finds the 
apparently quiescent phase likely composed of a signif- 
icant amount of microflaring emission, and by spectral 
analysis (Table [2]) which shows that active-region emis- 
sion persists throughout the duration of the observation 
at approximately the same temperature as found for the 
pre-flare quiescent emission (~ 5 MK). 

Together with the previously derived temperatures and 
emission measures, the observed flare rise time and the 
decay time scales can be used to determine the sizes and 
densities of the emitting plasma. When conduction losses 
dominate the fl are, the decay time scale can be appr oxi- 
mated (see e.g.. lWinebarger. Warren, fc Seatonll2003f ) by 

4xlO- 10 n e L 2 
Tc = ( CgS > ' W 

where n e is the electron density, L is the length scale over 
which the plasma is distributed (the flaring loop half- 
length), and T is the measured temperature. Given that 
the emission measure EM m \r? e l? for a hemispherical 
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volume of radius r = — , we can simultaneously solve 
for n e and r by identifying the observed decay time scale 
with tc ■ Similarly, when radiative losses dominate, the 
decay time scale is given by 

3kT 

TR = ^A(T)> (2) 

where A(T) is the intensity per emission measure (in erg 
cm 3 s" 1 ), obtained here combining line emissivities from 
ATOMDB (vl.3.1; Smith et al. 2001b) and continuum 
emissivities from SPEX (vl.10; Kaastra et al. 1996) with 
PINTofALE (Kashyap & Drake 2000) while using the 
fitted element abundances. As above, the plasma density 
n e and the radius of the emitting volume r can be derived 
from the measured temperature T and emission measure 
EM, after identifying the observed decay time scale with 

TR- 

At the temperatures typical of flares, A(T) oc T 1 ' 2 and 
so tr/tc oc T 3 . Conduction losses thus dominate during 
the early (hotter) part of the flare, while radiative cooling 
dominates the later phase. We compute estimates of the 
flare volume, plasma density, and loop sizes as follows: 

1. Flare Rise: The sound speed in a plasma at tem- 
perature T is c s = y 'kT j ' \xm v , where m p is the 
proton mass and \x is the average plasma particle 
mass coefficient (~ 1/2 in a fully ionized plasma). 
The rise time of a flare is generally due to the time 
it takes for evaporated chromospheric material to 
fill the flare volume. For the simplest case of impul- 
sive heating and a single loop, the rise time equals 
L/c s , and from this we can estimate the length 
scale of the flaring volume as L ~ 4 x 10 10 cm, or 
~ 3i?*, where i?* « 0.20 Rq is the radius of the 
star (cf. iSegransan et aTll2003[ ). Because L cx \/T, 
this result is only weakly sensitive to uncertain- 
ties in the flare temperature and provides a firm 
upper limit. In reality, the heating may be non- 
impulsive so that evaporation continues to fill the 
loop for longer than L/c s . The flare is also likely 
to occur in a complex loop arcade, in which the 
rise time is dominated by the successive lighting 
up of adjacent loop tubes (jReeves fe Warrenll2002l : 
iReeves et al.l 120071) . Solar loops have been mod- 
eled using hundreds of such strands, each lighting 
up a few seconds after the previous one; the typical 
velocities of the propagation of t he sideways distur- 
bance is ~ i the sound speed (jReeves fc Warrenl 
I2002D . With plausible assumptions regarding the 
geometry, this suggests typical flaring length scales 
of L ~ 4 x 10 9 cm. More accurate calculations 
require the use of the known decay timescales (see 
below) . 

2. Initial Decay: 

As noted above, temperatures in excess of 43 MK 
occur near the flare peak, and at such high tem- 
peratures, conductive loss is an important factor 
in the evolution of the flare intensity. Setting 
the observed decay timescale (r/ as t = 400 s) to 
match the conductive decay timescale (tc) in equa- 
tion 1 and using the observed flare emission mea- 
sure (EM w | 
sity n e ~ 1.6 x 10 12 



3 x 10 9 cm (corresponding flare volume of diameter 
~ 3.5 x 10 9 cm), in good agreement with the size 
derived from the flare rise time assuming an arcade 
of loops. Uncertainties on these results are large, 
however, because of the L cx T 5 /EM relationship 
and the large uncertainty (particularly on the high 
side) for the flare temperature. 

The expected radiative decay timescale for the de- 
rived plasma density is tr w 560 s, which is sim- 
ilar in magnitude to tq. Equating the two decay 
timescales (see van den Oord et al. 1988; their equa- 
tion 18), we estimate the loop semi-length as 



L«2.5x 10" 



EM 1 



cm , 



where Ni is the number of loops and a is the ra- 
tio of the loop diameter to its length, and where 
we have adjusted the coefficient by a factor of 0.7 
to account for the change in radiative power given 
the measured abundances. Assuming Ni = 1 and 
a = 0.1, we find L f=a 10 10 cm, greater than the esti- 
mate derived above. Modeling the flare decay more 
generally as a quasi-static cooling (van den Oord 
& Mewe 1989; their equation 15) results in a lower 
estimate of L ps 3 x 10 9 cm, matching the origi- 
nal simple estimate, again assuming a single flaring 
loop with a constant cross-section, and a = 0.1. 

3. Radiative Decay: As the flaring plasma cools, con- 
duction losses decrease in importance and radiative 
losses start to dominate (tr « tc). In the so- 
lar case, this phase coincides with the emergence 
of arcades of post-flare loops defined by the re- 
connected magnetic fields that cover approximately 
the same area as the pre-flare active region. Using 
equation (2), t r = t s1ow = 13000 s, T « 22 MK 
(kT = 1.9 keV), the observed decay-phase emis- 
sion measure (2.7 x 10 51 cm" 3 ), and tabulated val- 
ues for A(T), we derive a hemispherical volume 
of radius r ~ 7 x 10 9 cm filled with a plasma 
of density n e ~ 4 x 10 10 cm -3 . This is a ro- 
bust result, as L is fairly insensitive to temperature 
(L cx (EM/T) 1 / 3 ). For comparison, the pressure 
scale height during this phase, given by kT / (fj,m p g) 
and with the gravitational acceleration g calculated 
at the stellar surface assuming M* = O.2OM0, is 
~ 2.6 x 10 10 cm, and an arcade of magnetically 
confined loops that are one-quarter this height is 
quite plausible. 

Note that Reale, Serio, & Peres (1993; see also 
Reale et al. 2004) have developed a modeling 
method that fits hydrodynamically evolving loops 
to the observed light curves. Based on simulations, 
they derive a scaling law to determine the loop 
semi-length based on the initial temperature of the 
flare and the decay time scale, 



1 n 2 L 3 ), we derive a plasma den- 



cm and loop size L 



x/T/lOpVIKl o 

Lwr *"» 120/(0 10 

where /(£) > 1 is a non-dimensional correction 
factor that takes into account whether heating is 
present during the decay and £ is the slope of 
the decay path in a density-temperature diagram. 
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Note that C ~ 1-2 here (see jLl below). Us- 
ing this scaling law, we derive a loop semi-length 
L ~ 7 x 10 s cm during the fast decay. Ex- 
trapolating to the slow decay phase, we derive 
L ~ 2 x ^^ 10 7TC) cm ' ^ e ^ a ** er * s consistent 
with the estimates derived above, but the former 
greatly underestimates the length scales involved. 
Not e that modeling this decay as quasi-static cool- 
ing (jvan den Oord fc Mewdll989t) results m an es- 
timate of the loop semi-length L w 3 x 10 11 cm, 
which is greater than all other estimates as well as 
greater than the coronal pressure scale height. We 
thus conclude that the slow decay is not well ex- 
plained by quasi-static cooling, and that the simple 
picture of a single hydrodynamically evolving loop 
is too simplistic to explain the characteristics of 
this flare. A detailed hydrodynamic modeling of 
this flare is in progress. 

4. Active Region: In quiescence (i.e., prior to the 
flare) , the active region(s) on the star cover an area 
w 50 ± 13% (Johns-Krull & Valcnti 1996) of the 
stellar surface. Despite the dynamic nature of an 
active region, based on solar observations, we may 
approximate it as a stationary steady state loop in 
radiative equilibrium (e.g., Rosner, Tucker, & Va- 
iana 1978 [RTV]). In such a case, the temperature 
T, pressure p ps 2n e kT, and the loop length L obey 
the scaling law 

T = K (pL)i, (3) 

where the proportionality constant is estimated to 
be k ~ 1.4 x 10 3 in the case of the Sun. The value 
of k is dependent on the assumed estimate of the 
radiative loss function. For the case of Ross 154, we 
estimate that k is reduced by a factor of 2.05±0.20. 
The emission measure of such a plasma is given by 

EM = n 2 e V w n 2 e (fAnRlL) , (4) 

where / is the fraction of the stellar surface covered 
by active regions. For / = 0.5, and the measured 
emission measure of the quiescent emission at T = 
5.3 MK (kT = 0.46 keV), we find loop lengths L ks 
4 x 10 10 cm. This is much larger than both the 
loop sizes estimated above and the pressure scale 
height at this temperature (~ 6 x 10 9 cm). The 
latter violates the RTV assumption of a ^constant 
pressure environment, and therefore we conclude 
that the active region loops are probably not in 
a stationary state. L, however, is proportional to 
t^T JEM , and at the assumed temperature K is 
approximately proportional to T~ 2 so that L cx 
T~ 8 /EM. A 10% uncertainty in T therefore leads 
to an uncertainty of ~ 2.5 in L. 

The size and density estimates derived above are all 
subject to significant uncertainties, because of assump- 
tions regarding geometry (arcade versus single loop for 
the Flare Rise, low-lying stationary loops for the Active 
Region) or uncertainties in the fitted temperatures (par- 
ticulary in the Initial Decay analysis, but also for the Ac- 
tive Region and Radiative Decay calculations). Results 
from the various flare analyses are reasonably consistent 
but we hesitate to interpret them in more detail. 



7.2. Color-Intensity Evolution 

Detailed spectral analysis (e.g., Tabled]) can only be 
carried out over durations long compared to the time 
scale over which the flare evolves; analyses over shorter 
durations result in large errors on the fit parameters be- 
cause of inadequate counting statistics. We therefore 
track the evolution of the hardness ratio of the flare as a 
means to study the evolution of the flare plasma. This 
evolution is shown in a color- intensity plot (Figure 1 1 0[) , 
where the color C is the log of the ratio of counts in 
the soft (S; 0.25 - 1 keV) and hard (H; 2 - 8 keV) pass- 
bands and the intensity r is the combined count rate 
in those two bands. We compute the hardness ratios 
and their errors using the full Poisson likelihood in a 
Bayesian context (Park et al. 2006). This method allows 
us to determine the error bars accurately even in the low- 
count limit, and to explicitly account for and incorporate 
the systematic uncertainty inherent in a specific choice 
of time-bin size St and time-start phase by considering 
all possible choices of phase (via cycle-spinning) and bin 
sizes (via marginalization). The resulting track of the 
flare is shown in Figure I10[ where darker shades repre- 
sent a longer time spent by the source in that part of 
the color-intensity space and the fuzziness of the shading 
reflects the statistical error in the calculation. 

Also shown on the plot arc a set of line segments 
marking the evolutionary track of the flare for the spe- 
cific time-bin size St = 150 s, and a pair of model 
grids computed for different values of emission mea- 
sure that bracket the observed flare intensities (EM = 
0.5 — 15 x 10 52 cm" 3 ). The grids are computed with 
the PINTofALE package (Kashyap & Drake 2000) us- 
ing atomic emissivities obtained with the CHIANTI v4.2 
package (Young et al. 2003) and ion balance calculations 
by Mazzotta et al. (1998). 

A number of notable features are present in the evolu- 
tionary track of the flare. First, there is a sharp initial 
hardening of the spectrum coincident with the onset of 
the small initial optical flare (bins 1-3; see Figures [9] 
and ITU|) . We identify this with the initial onset of the 
reconnection event that generates a non-thermal hard 
X-ray flare. The spectrum at bin 3 appears to be ei- 
ther non-thermal so that the EM grids do not apply, or 
else has a high temperature and a very high column den- 
sity (Njj > 5 x 10 20 cm" 2 ) that might be associated 
with a coronal mass ejection. There are, however, too 
few counts for a spectral-fitting analysis to distinguish 
between these possibilities, and the uncertainty ranges 
for r and C may be large enough to permit a less in- 
teresting low-density thermal explanation. Thereafter, 
the plasma thermalizes, causing the spectrum to become 
softer and the intensity to rise because of increased radia- 
tive loss (bins 3-8). At this point, the main reconnection 
event occurs, leading to a large flare that increases in 
intensity and spectral hardness (bins 8-11), before heat 
input ceases and the plasma decays back to lower inten- 
sities and softer spectra (bins 12-19) during the rapid 
conductive-decay phase of the flare. 

The color- intensity track in Figure HU1 can be converted 
to a locus of evolution in Temperature and Emission 
Measure (T — EM) space if the absorption column Nh 
and the plasma abundances are fixed throughout the du- 
ration of the flare. We adopt a low column density of 
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Fig. 10. — Color-intensity evolution of the flare, obtained by 
averaging cycle-spun images at various bin sizes ranging from 40 s 
to 400 s. The density plot is made by sampling from the posterior 
probability distribution of the color C (the log of the hardness ratio 
S/H) computed from the binned data points for various bin sizes 
and is also cycle spun to include the effects of varying the phase 
of the binning. Data points obtained for a light curve constructed 
for a bin size of 150 s (see Fig. [9} are marked on the plot and 
connected by the solid white lines to trace the evolution of the 
source. A clear pattern in the flaring plasma evolution can be 
seen, first rising in temperature and intensity, softening at the peak, 
and then decaying in temperature and intensity. Also shown are 
grids corresponding to an optically-thin thermal emission model 
calculated for two different emission measures, 1.5 X 10 53 cm -3 for 
the upper one, and 5 X 10 51 cm -3 for the lower one. The grids 
are computed for temperatures T = 10, 20, 50, and 90 MK and for 
column densities log ATjy = 19 to 21.3. Spectral fits to the summed 
data result in values of Njj that are consistent with zero, however, 
from comparing the model thermal emission grids with the color- 
intensity values for bin 3, instances of large Nh cannot be ruled 
out. 
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Fig. 11. — The evolutionary track of the main flare in 
Temperature-Emission-Measure space. A plasma model assum- 
ing the same abundances as determined with the spectral analysis 
of the flare counts (Table |2j and a low column density (Njj = 
10 1 9 cm " 2 ) was used to convert the colors and intensities of Fig- 
ure 1101 to temperature and emission measure. The bin numbers 
corresponding to the time of the flare (see Figure |9j are shown, 
connected by dotted lines. 

Njj = 10 19 cm -2 (any value below ~ 10 20 cm -2 has no 
significant effect) and the fitted flare abundances (Ta- 
ble [2|) to compute the temperature from the color and 
the emission measure from the intensity (Figure [TTj) for 
the main flare (bins 8-19). If we further assume that 
the flare volume does not change over this time, then the 
EM axis tracks changes in n 2 . From the T — EM track, 
we see this sequence of events: first, the temperature in- 
creases rapidly, followed by a rapid increase in density, as 



the evaporated chromospheric material fills the flare vol- 
ume (bins 8-11). Then, the plasma starts to cool rapidly 
as the density reaches its peak, and eventually returns to 
the pre-main- flare environment. 

During this cooling phase (bins 12-19), EM oc 
T 2 4±0 - 7 . Assuming that the flare volume does not vary, 
this means that n e oc where £ w 1.2 ± 0.35 (see Ra- 
diative Decay discussion in $7.1\i . This value of ( corre- 
sponds to an intermediate level of heating occuring dur- 
ing the flare, similar to the flare observed on Proxima Cen 
(Reale et al. 2004). The flare evolution, however, is seen 
to be highly complex; looking at a shorter portion of the 
cooling phase (bins 13-16), the index £ « 0.6±0.1, which 
is indicative of greater heat deposition. This continuous 
heating may be a manifestation of a flare arcade, as is 
often seen during solar flares (e.g., Reeves et al. 2002). 
Detailed hydrodynamical modeling of the flare arcade is 
necessary in order to definitively establish the flare envi- 
ronment. 

8. MASS LOSS AND STELLAR WIND CHARGE EXCHANGE 

HALO 

Mass-loss rates of a few times 10~ 10 to more than 10 -5 
Mq yr _1 have been measured for many stars using a va- 
riety of methods based on P Cygni profiles, optical and 
molecular emission lines or absorption lines, a nd infrared 
and radio excesses (|Lamers fc Cassinc lli 1999). All those 
measurements, however, are of strong winds from mas- 
sive OB and Wolf-Rayet stars or cool red giants and su- 
pergiants. For comparison, the solar mass-loss rate {Mq) 
is only ~ 2 x 10 -14 Mq yr -1 . Because of their frequent 
large flares, it might be expected that dMe stars would 
have comparable or larger mass-loss rates, and because 
of their sheer numbers, M dwarfs would then contribute 
substantially to the chemical enrichment of the ISM. 

It was not until very recently that measurements of the 
modest stellar winds from low-mass stars became feasi- 
ble, beginning with a Cen AB (M= 2M Q ) and Proxima 
Cen (M< O.2M ) (jWood et al.ll2001h . Over a dozen late- 
type mostly main-sequence stars have now been stud- 
ied with measured mass-loss rates ranging from 0.15 to 
100 (iMuller. Zank. fc Woodll200l iWood et al.ll2002l . 
|2005|) . Those studies have revealed that mass- loss rates 
arc roughly proportional to a star's X-ray luminosity per 
unit surface area, except for M dwarfs, which surpris- 
ingly have much lower mass-loss rates than predicted by 
that relation. 

All the late-type-star measurements were made using 
the Hubble Space Telescope Goddard High Resolution 
Spectrometer to measure H Lya absorption profiles. To 
summarize, the interaction of a stellar wind with a par- 
tially neutral ISM, primarily in charge exchange (CX) 
collisions between stellar-wind protons and neutral H 
atoms from the ISM, creates a region of enhanced neu- 
tral H density, or 'hydrogen wall,' in the star's astropausc 
(analog of the Sun's heliopause). The resulting density 
enhancement leads to excess absorption in the wings of 
the star's Lya line, which can be measured and modeled 
to deduce the total stellar mass-loss rate. 

At the time the first of those results were published, 
another more direct but less sensitive stellar-wind detec- 
tion method, al so based on charge excha nge processes, 
was proposed bv lWargelin &: Drake! (|2001h . This method 
searches for X-ray emission from the very small fraction 
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of metal ions in the stellar wind, which charge exchange 
with neutral gas streaming into the astrosphere from the 
ISM. In CX collisions of the highly charged wind ions, 
especially H-like and fully ionized O, an electron from a 
neutral H or He atom is captured into a high-n level of 
the metal ion (usually n = 5 for O) from which it then 
radiativcly decays, emitting an X-ray. Although coro- 
nal X-ray emission from a star is roughly 10 4 times as 
bright as its stellar-wind CX emission, the latter is emit- 
ted throughout the astrosphere with a distinctive spec- 
tral signature which permits both spatial and spectral 
filtering to be applied. This X-ray CX method was first 
applied to a Chandra observation of Proxima Cen, result- 
ing in a null detection of the CX halo and an upper limit 
for the mass- loss rate of 14M () Wargelin fc Drakdl2002f ). 
compared to an upper limit of 0.2M W derived using the 
Lya method ([Wood et al.ll200ll ). In the remainder of §8 
we apply the CX technique to Ross 154, explain why we 
are ultimately unable to deduce even an upper limit for 
its mass-loss rate, and discuss future prospects for this 
detection method. 

8.1. CX Analysis 

As described by IWargelin fc Drake! (|2002ft , the total 
rate of X-ray emission from a star as the result of CX 
of a particular species of ion is simply equal to the pro- 
duction rate of that ion in the stellar wind, either from 
its creation in the corona or as the result of the CX of 
a more highly charged ion state. As an example, the 
emission rate of He-like O X-rays is equal to the sum 
of the creation rate of H-like O ions and the creation 
rate of bare O ions, which eventually charge exchange to 
create H-like ions. Oxygen is the most abundant metal 
ion in stellar winds, and with a coronal temperature of 
0.46 keV (5.3 x 10 6 K) nearly all of the O ions in Ros s 
154's wind are fully stripped 8+ (jMazzotta et al.lll998l ). 
Those ions emit O VIII Lyman photons via CX, followed 
eventually by O VII Ka emission after the ions CX a 
second time. 

The total number of X-ray photons emitted via CX 
processes from oxygen ions in the Ross 154 wind and 
then detected by a telescope with effective area A is thus 
equal to 

where Ro is the production rate of (fully stripped) O 
ions in the stellar wind, t is the observation time (42620 
s for the quiescent phase), and d is the distance to the 
star (2.97 pc). The factor of 2 is because each O ion 
emits a H-like and then a He-like X-ray. 

To calculate Ro we consider the case in which Ross 154 
has the same mass-loss rate as the Sun, M= 2 x 10~~ 14 
Mq yr _1 . Assuming the same He/H ratio as for the so- 
lar wind, 0.044, and half the metal abundance (net 0.001 
relative to H), the total mass-loss rate is 6.65 x 10 35 ions 
s _1 , of which 6.36 x 10 35 ions s _1 are protons. The ratio 
of O and H ions in the wind (which we assume is equal to 
their ratio in the corona) is equal to the product of: the 
coronal abundance of O relative to the solar photosphcric 
abundance; the s olar photospheric to coro nal O abun- 
dance ratio (T.10; IAnders fc Grevessell 1989ft: and O to H 
ratio in the solar wind f5.6x 10~ 4 : ISchwadron fc Cravens! 
[2000]) . In the case of half-solar coronal abundances, Rq is 



1.9 x 10 32 ions s , and with an effective area of ^320 cm 2 
for C7ian(ira/ACIS-S3 in the energy range for O emission 

(~600 eV), N t ot is 5 counts for a 1M Q wind. If those 
counts could be isolated from the detector background 
and coronal emission, a stellar wind only a few times 
stronger could be detected. 

8.1.1. Searching for the CX Excess 

To isolate the CX emission we take advantage of its 
unique spectral and spatial characteristics. As noted 
above, oxygen is the most important element in stel- 
lar CX halos, primarily because it is the most abundant 
metal. It is also the highest-Z element that has a large 
fraction of the bare and H-like ions that emit X-rays fol- 
lowing CX; C ions are nearly abundant and even more 
highly ionized, but their X-ray emission is at lower en- 
ergies where X-ray detection efficiency is usually poor. 
L-shell Fe emission may also contribute but at a lower 
level over a relatively broad range of energies. For these 
reasons we focus on O emission between 560 and 850 
eV, encompassing He-like O VII Ka (^565 eV), H-like 
O VIII Lya (654 eV), and the high-n Lyman lines (775- 
836 eV). Because of the moderate energy resolution of 
the ACIS-S3 CCD, we make our energy cuts at 525 and 
875 eV. 

To further enhance contrast against the coronal emis- 
sion we use spatial discrimination and compare the radial 
profiles of the source emission during quiescence and dur- 
ing the flare (when the fraction of CX emission is negli- 
gible) . As mentioned above and explained in more detail 
in £18.1.21 Ross 154's CX emission is expected to appear 
around the star in an approximately symmetric shell- 
like halo. During the flare the observed source emission 
distribution will be essentially identical to the intrinsic 
telescope PSF because any CX emission will contribute 
a minuscule fraction of the total observed events. Dur- 
ing quiescence any halo emission will, if strong enough, 
manifest itself as an excess relative to the the flare PSF. 
Although the flare and quiescent spectra are different, 
the energy range we use for our analysis (525-875 eV) 
is narrow enough that differences in the effective PSF 
for each phase's emission are negligible compared to the 
statistical uncertainties described below. 

When searching for the CX halo emission, we cannot 
look too close to the source centroid because the coronal 
emission is thousands of times stronger and completely 
swamps that from the halo, despite the very tight Chan- 
dra PSF. Likewise, the halo extraction region cannot be 
too large or the detector background will dominate the 
CX signal. For this observation (with its particular back- 
ground, and number of quiescent and flare counts within 
the chosen energy band), an annulus with radii of ap- 
proximately 35 and 81 pixels (corresponding to 17" and 
40", or 51 and 119 AU) provides the best sensitivity af- 
ter a 5-pixel-widc box around the bright readout streak 
is excluded. 

To compare the fraction of halo-region emission during 
the quiescence and flare we must know the total number 
of counts during each phase. This was determined from 
the number of counts in the readout streaks (excluding 
a 20-pixel-radius circle around the source and correcting 
for background and underlying events distributed by the 
telescope PSF), dividing by the fraction of the streak 
used (162 rows out of 206), and then dividing by the 
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exposure time fraction for the complete readout streak, 
equal to 1.35% (see p~2J|) . 

Using the energy range and source region described 
above we found 97 ± 14 quiet counts and 200 ±15 flare 
counts, giving count fractions of 0.00285 ± 0.00044 and 
0.00367 ± 0.00032, respectively. The net difference is 
therefore -0.00082 ± 0.00054. Equivalently, 125 ± 13 
quiet counts were expected based on the flare ratio while 
97 ± 14 counts were observed, yielding a statistically in- 
significant deficit of 28 ± 19 counts. Several other combi- 
nations of energy range and annulus size were also tried, 
with similar results. 

8.1.2. Modeling the Spatial Distribution 

After obtaining the above null result for CX halo emis- 
sion, we attempted to obtain an upper limit on the stellar 
mass-loss rate by computing how many CX events would 
need to be detected to obtain a statistically significant 
excess. This requires more detailed modeling of the ex- 
pected spatial distribution of the CX emission in order to 
determine / onn , the fraction of CX photons that would 
be detected within the chosen annulus. 

Based on absorption-line spectra obtained by HST and 
EUVE, Ross 154 is believed to be surrounded by the 
'G cloud,' a r egion of enhanced de nsity a few parsecs in 
size for which iLinskv et alJ (|2000l ) have derived a value 
of 0.10 cm~ 3 for the neutral H density. The G cloud 
is moving roughly toward us from t he direction of Ross 
154 toward I = 184.5°, b = -20.5° (jLallement fc Bertinl 
[l992f h with virtually the same velocity (29.4 km s _1 ) 
and direction of motion as the Local Interstellar Cloud 
surrounding the Sun. Given its proper motion and the 
radi al velocity listed in the SI MBAD database (-4 km 
s" 1 : IWallerstein fc Tvagill2004h . Ross 154 sees an ISM 
wind of 22 km s _1 , with the upwind side of its astro- 
sphere pointed roughly toward the Galactic center and 
the viewing angle from Earth nearly aligned with the axis 
of its astrosphere (at an angle of 26°; B. E. Wood 2005, 
private communication) . 

As noted in SJ5J a hydrogen wall of neutral gas is ex- 
pected to form on the upwind side, and it is in this region 
of enhanced neutral gas density that the CX emission 
will be most concentrated. Given the orientation of the 
H wall, we model it as a hemispherical shell pointed to- 
ward Earth. Outside the shell, the ISM is undisturbed, 
while inside the density of neutral H can be approximated 
(|Cravensll2000l ) as nn 



•h 



-(•WO where 



asy • 



is the 



asymptotic density just inside the H wall, A# is the scale 
for neutral H depletion near the star (from photoioniza- 
tion and proton-H CX), and r is radial distance. 

Based on the astrospheric modeling results of 
iWood et all (p002Tl we have derived the following approx- 
imate scaling relations: 
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where Rmner and R ou ter are the inner and outer radii 
of the H wall, n™ n is the neutral H density within the 
wall, njf M is the neutral H density in the undisturbed 
ISM (0.10 cm~ 3 ), and all values are appropriate for the 
upwind direction for M/Mq> 1 and relative cloud/star 
velocities of less than ~ 50 km s _1 . 

The distribution of neutral He is less well understood 
but also less important. We assume that its density is 
given by iijj e = O.lOn^ ' e -( x H<;/ r ) where A# e is approx- 
imately 1 AU. CX emission from the downwind 

hemisphere is roughly equal to the upwind emission but 
spread over such a large volume of space that we ignore 
it in our detectability calculations. 

Detailed descriptions of the spatial emiss io n cal- 
culations arc found in IWargelin fc Drake] (|2001l . 
l2002f >. Cross sections for CT of highly charged 
C, N, O, and Ne ions with atomic H and asso- 
ciated radiative br anchi ng r atios are taken from 
lHarel. Jouin. fc Ponsl (|1 998ft. [Greenwood et all 
2001h. iRigazio. Kh archcnko. fc Dalgarnol (|2002f ). 

„i do 



Johnson et ajj (I2000D. and related references m 
Kharchenko fc Dalgarnol (|2000L 1200 If ). To summa- 



rize, the stellar wind is assumed to expand isotropically 
and CX with neutral gas in the astrosphere having a 
distribution and density described by the equations 
above. For bare 8+ and H, the CX cross section acx 
is 3.4 x 10~ 15 cm 2 , and 3.7 x 10" 15 cm 2 for 7+ . CX 
with He, which is relatively unaffected by its passage 
from the ISM into the astrosphere (there is no 'He wall' 
density enhancement) is also included in our model. 
The density of each ion species is calculated numerically 
and its emissivity and fractional abundance computed 
at each step (e.g., the fraction of 8+ decreases with 
distance as it charge exchanges to create 7+ ions and 
O VIII Lyman photons). The initial composition of 
stellar- wind O ions, is set to 95% 8+ and 5% 7+ . 

For modest M / M@ (less than ~20), the CX emission 
profile (projected on the sky and summed in annular bins 
around the star) rises steeply from the center and then 
slowly falls (sec Fig. [12]). The deficit of emission at small 
radii is because of neutral gas depletion near the star, 
and the falloff in X-ray CX emission at large radii occurs 
because highly charged ions arc 'used up' as they CX 
with neutral gas, recombining into lower charge states 
that can not emit at X-ray energies. The CX luminosity 
from such winds, however, is too small for us to detect. 

For larger mass-loss rates, the astrosphere is essentially 
swept clear of neutral gas {n^ v < 0.10; see eq. 9) so that 
CX emission is concentrated in the astropausc in a shell 
around the star. When viewed from afar the emission 
per radial bin steadily rises from the center as the line 
of sight traverses an increasingly long path through the 
emission shell and then gradually falls at larger distances 
as the fraction of highly charged ions decreases. The 
emission shell is hundreds or thousands of AU in radius 
in such a case (eqs. 6 and 7), while the annulus used in 
our analysis extends to only 119 AU because of the lim- 
ited field of the observation (and because the background 
per unit radius becomes too large). The fraction of total 
CX emission within this annulus, / arm , therefore scales 
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Fig. 12. — Modeled radial CX emission distributions for various 
mass-loss rates. 100 AU corresponds to 34". At rates above ~ 
20Mq, most emission occurs in a shell hundreds of AU from the 
star, and the emission observed in an annulus near the star (51—119 
AU, shaded gray) becomes increasingly independent of the mass- 
loss rate. Emission at small radii is nearly all from H-like 7+ but 
at larger distances an increasing fraction of emission comes from 
Hc-likc O e+ as stellar-wind ions undergo further CX collisions. 

as approximately Runner- According to equation 6, how- 
ever, Rinner is proportional to A/ 5 , so f ann scales as 
M -1 . Since the total number of counts detected within 
the annulus is given by R 

ann — farinRtot? and Rtot is 

proportional to M, R an n is nearly constant. Our ability 
to detect the CX excess therefore does not increase with 
M beyond ~2OM , and so we are unable to find even an 
upper limit. In short, a CX halo provides too few counts 
to be detectable for small M, while larger values of M 
produce more total emission but spread too diffusely over 
a larger projected area. 

8.1.3. Future Prospects for CX Halo Detections 

Ultimately, the limiting factor in such a search is tele- 
scope sensitivity. A large collecting area is thus essen- 



tial, while detectors with better energy resolution can 
vastly improve the contrast of the CX signal against 
coronal emission and background. Micr ocalorimeters, 
such as the now-lost XRS onboard Suzaku (|Cottam et al.l 
l2005h and those planned for proposed missi ons such as 
Constellation- X (IWhite fc Tananbaumll2003f ) and XEUS 
(|Parmar et al.ll2003[ ). have resolution of a few eV and can 
isolate individual lines that are enhanced by CX, such 
as high-n Lyman lines and the He-like Kg forbidden and 
interc ombination lines jWargelin . Beiersdorfer. fc Brownl 
12001 . With collecting areas roughly 10 and 100 times 
that of Chandra, Con-X and XEUS could respectively 
collect hundreds or thousands of CX counts from solar- 
mass-loss-rate stars in 100-ks exposures. 

Such observations would permit direct imaging of the 
winds and astrospheres of main sequence late-type stars 
(i.e., stars with highly ionized winds), providing informa- 
tion on the geometry of the stellar wind, such as whether 
outflows are primarily polar, azimuthal, or isotropic, and 
whether or not other stars have analogs of the slow (more 
ionized) and fast (less ionized) solar wind streams. En- 
hanced emission in the relatively dense neutral gas out- 
side the stellar-wind termination shock would indicate 
the size of the astrosphere and the orientation and ap- 
proximate speed of the relative star-cloud motion. Neu- 
tral gas density can be inferred from the falloff of CX 
emission at large distances from the star, which is pro- 
portional to 1/tcx = 1/ '^cx^h, where 1/rcx is the 
path length for a CX collision. The stellar wind velocity 
can also be determin ed from the hardness rati os of the H- 
like Lyman emission ([Beiersdorfer et aLll2001[ ) . Although 
such studies are just beyond the capabilities of Chandra, 
future observatories will certainly provide abundant in- 
formation on the stellar winds and astrospheres around 
nearby stars. 

9. SUMMARY 

In this paper we have analyzed data from two Chandra 
observations of Ross 154, which included one very large 
and a few moderate flares. Analysis of spectra from the 
ACIS observation reveals an enhanced Ne/O abundance 
ratio (relative to commonly accepted solar abundances) 
during quiescence and a further increase, along with over- 
all metallicity, during the large flare. Because of severe 
event pileup, however, large corrections to the instrument 
effective area were needed and each spectral fit used only 
a few thousand counts, introducing significant uncertain- 
ties into the results. 

Flaring behavior was studied at both low and very 
high intensity levels. Two moderate flares, which sat- 
urated the limited available telemetry, occurred during 
the HRC-I observation, and evidence for the Ncupcrt ef- 
fect was seen at the beginning of the large ACIS flare. 
Based on the spectral fit results and measured flare rise 
and decay times, we derived several physically reason- 
able estimates of flare and active region sizes and den- 
sities, along with information about the color-intensity 
evolution of the flare. 

From temporal analysis of the quiescent emission dur- 
ing both observations we found that microflaring, which 
is the superposition of many low-level flares, shows a 
bimodal intensity-frequency distribution when modeled 
with a power-law distribution. Only a few other stars 
have been similarly studied, and our results provide the 
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first evidence that flaring on late-type stars may follow a 
distribution more complex than the single power-law of 
the Sun. Lastly, we searched unsuccessfully for charge- 
exchange emission from the stellar wind, but conclude 
that the next generation of X-ray telescopes will be ca- 
pable of such observations. 
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